KEK Preprint 2010-54 
IPMU11-0012 
TU-875 
UT-11-18 

Discrimination of New Physics Models 
with the International Linear Collider 

Masaki Asano^ a )0, Tomoyuki Saito^, Taikan Suehara( c fl 
Keisuke Fujii^, R. S. Hundi^, Hideo Itoh^, Shigeki Matsumoto^, 
Nobuchika Okada^, Yosuke Takubo^, and Hitoshi Yamamoto' 6 ' 

( a ) Department of Physics, The University of Tokyo, Tokyo 113-0033, Japan 
^ Department of Physics, Tohoku University, Sendai 980-8578, Japan 
^ICEPP, The University of Tokyo, Tokyo 113-0033, Japan 
^High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801, 

Japan 

^ Department of Theoretical Physics, 
Indian Association for the Cultivation of Science, Kolkata 700032, India 

^ICRR, The University of Tokyo, Kashiwa 277-8582, Japan 
^IPMU, TODIAS, The University of Tokyo, Kashiwa 277-8582, Japan 
(h) Department of Physics and Astronomy, University of Alabama, 
Tuscaloosa, AL 35487, USA 



1 masano@hep-th.phys.s.u-tokyo. ac.jp 
2 saito@epx. phys.tohoku.ac.jp 
3 suehara@icepp.s.u-tokyo. ac.jp 



Abstract 

The large hadron collider (LHC) is anticipated to provide signals of new 
physics at the TeV scale, which are likely to involve production of a WIMP 
dark matter candidate. The international linear collider (ILC) is to sort out 
these signals and lead us to some viable model of the new physics at the TeV 
scale. In this article, we discuss how the ILC can discriminate new physics 
models, taking the following three examples: the inert Higgs model, the super- 
symmetric model, and the littlest Higgs model with T-parity. These models 
predict dark matter particles with different spins, 0, 1/2, and 1, respectively, 
and hence comprise representative scenarios. Specifically, we focus on the pair 
production process, e + e" — > X + X~ ~^ X°X°W + W~ , where x° an d X^ are 
the WIMP dark matter and a new charged particle predicted in each of these 
models. We then evaluate how accurately the properties of these new parti- 
cles can be determined at the ILC and demonstrate that the ILC is capable 
of identifying the spin of the new charged particle and discriminating these 
models. 
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1 Introduction 



The standard model (SM) has no symmetry to protect the smallness of the scale of 
the electroweak symmetry breaking, hence, the Higgs mass receives quadratically di- 
vergent corrections, leading to the hierarchy problem. As a remedy for this problem, 
new physics beyond the SM is expected to appear at the TeV scale. 

The SM has, however, yet another problem. We know that about 23% of the 
energy density of the present universe is made up of unknown dark matter [1] and 
that it played an important role in the formation of the large scale structure of the 
universe [2]. There is, however, no candidate for the dark matter in the SM. 

It seems plausible that the problem of the dark matter is also solved in the 
framework of the physics beyond the SM which solves the hierarchy problem. In the 
TeV scale physics, there are new particles which change the behavior of the quantum 
correction to the Higgs mass term. Some of the new particles would have the SU(2) 
charge of the SM because it is related with the origin of the electroweak symmetry 
breaking. To solve the hierarchy problem without a fine-tuning, these have masses 
of (9(100) GeV. When the lightest of them is neutral and stable (e.g. the lightest 
neutralino in supersymmetric models with conserved R-parity), it is nothing but a 
Weakly Interacting Massive Particle (WIMP) [3]. The WIMP is well known to be a 
good candidate for the dark matter, which naturally realizes the correct dark matter 
abundance in the present universe. Because of these attractive features, many new 
physics models at the TeV scale involving the WIMP dark matter candidate have 
been proposed. 

One of the most important questions here is how to single out the new physics 
model at the TeV scale that consistently describes the results from energy frontier 
colliders such as the Large Hadron Collider (LHC) and the International Linear 
Collider (ILC). Uncovering the nature of the WIMP dark matter experimentally is 
of particular importance not only for particle physics but also for astrophysics and 
cosmology. The LHC experiments are now in operation where new physics signals 
are anticipated, which will guide us to narrow down possible models at the TeV 
scale. Being a hadron collider, the LHC is suitable to study colored new particles. 
It is, however, not an ideal place to do precision measurements of the properties of 
weakly interacting particles (non-colored particles) including the dark matter, while 
the ILC, being a lepton collider, has a great advantage for this purpose. 

In this article, we investigate the possibility to discriminate new physics models 



1 



at the ILC using the following process including the WIMP dark matter: e + e _ — > 
X + X~ ~~ X°X°W + W~ , where x° an d X ± stand for the WIMP dark matter and a 
charged new particle predicted in each model. In our previous work [3], we have 
investigated the process in the framework of the littlest Higgs model with T-parity, 
and evaluated how accurately we can measure the properties of the new particles. We 
have shown there that the masses of X s1 an d X° can be determined to an accuracy 
of 1% or better by locating the both endpoints of the energy distribution of the 
reconstructed W bosons. It is also possible to determine the spin of x ± an d the 
structure of the interaction vertex between x ± i X°; an d W, through the observations 
of the angular distribution of x ± an d the polarization of W. The gauge charge of x^ 
can also be measured, making use of polarized electron beam. Interestingly, the same 
process exists in various other new physics models at the TeV scale, and it turns out 
to be an extremely useful process to extract information on the new physics. 

As the first step of our study to evaluate the ILC's potential to single out a viable 
new physics model, we investigate the possibilities to discriminate the following three 
models: the inert Higgs doublet model [5], the supersymmetric model [6], and the 
littlest Higgs model with T-parity |7j. These models contain a WIMP dark matter 
particle with spin 0, 1/2, and 1, respectively [8]. The masses of x an d X° are 
adjusted to coincide among different models. Although these models predict different 
cross sections for the x P arr production, we also force the cross sections to be a 
common value. We thus concentrate on the information related to the spin of the 
new charged particle for the discrimination of the new physics models. 

This article is organized as follows. In the next section, we briefly review the new 
physics models used in our simulation study. Simulation framework such as repre- 
sentative points and simulation tools is presented in section 3. Details of the analysis 
to discriminate the new physics models are given in section 4, where the expected 
measurement accuracies of x ± an d X° properties are shown for each representative 
point. Section [5] is devoted to summary. 

2 New Physics Models 

As already mentioned in the previous section, we concentrate on the process, 

e + e~ -> X + X~ -> X°X°W + W-, (1) 
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Particles 



Spins 



Representative Model 



(xh Xs) 
(x% X° F ) 
(xj, X°v) 
(Xy> Xs) 
(Xs, Xv) 



(0, 0) 
(1/2, 1/2) 

(1, 1) 
(1, 0) 
(0, 1) 



Inert Higgs model 
Supersymmetric model 
Littlest Higgs model 
No well-known models 



No well-known models 



Table 1: Spins of new particles x° an d X m various new physics models. 

where the WIMP dark matter is denoted by x°> while the new charged particle is 
X , and both particles are assumed to have odd charge under the Z 2 symmetry 
guaranteeing the stability of the dark matter. SM particles are assumed to have 
even charge under the symmetry. The interaction vertex between x°, x^i and W T 
exists in the most of new physics models at the TeV scale. On the other hand, the 
spins of x° an d X ± are dependent on the model. All possible combinations for the 
spins up to spin 1 are shown in Table [TJ 

At the ILC, x + an d x~ are produced in pairs through s-channel exchanges of 
photon (7) and Z boson, and the produced x ± decays to x° and W ± . In addition, 
if there is another new particle which has a lepton number such as the sneutrino in 
the supersymmetric model or the heavy neutrino in the littlest Higgs model with 
T-parity, the diagram in which the new particle is exchanged in the t-channel con- 
tributes to the x ± pair production. In our analysis, we simply assume that such a 
particle is heavy enough and ignore its contribution. 

In our simulation study, we consider the inert Higgs doublet model, the super- 
symmetric model, and the littlest Higgs model with T-parity as benchmark models 
in which the x ± has spin 0, 1/2, and 1, respectively, and develop the strategy to 
discriminate these models at the ILC. The crucial difference from the (1,0) or (0,1) 
models in Table [1] only appears in what relates to the x ± X°W T vertex (e.g. the 
shape of the energy distribution of W bosons), so that the strategy developed in this 
article can be applied to the models with (1,0) or (0,1) spin combinations. 

In the rest of this section, we briefly introduce the models used in our simulation 
study, focusing on interactions relevant to our analysis. 
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2.1 Inert Higgs doublet model 

The inert Higgs doublet model [5] is one of the two-Higgs-doublet models with un- 
broken Z2 symmetry. One of the Higgs doublets transforms as <f> •<->■ — under the 
discrete symmetry, while the other doublet and SM particles transform as SM -h- SM. 
Because of the existence of the terms which break the custodial symmetry in the 
Higgs potential, the mass of the lightest Higgs boson could be as large as 500 GeV 
without conflicting with precision electroweak measurements. The fine-tuning be- 
tween the Higgs boson mass and its radiative corrections, therefore, becomes mild 
compared to the SM with a light Higgs boson. In the model, neutral and charged 
components of the Z2 odd Higgs boson, which is called the inert Higgs boson, play 
the role of the WIMP (x°s) and the new charged particle (xs)> both of which are 
scalar particles. 

We focus on production and decay vertices of the new charged particle, which 
originate from gauge interactions: 

C = 1 [g z (l/2 - s 2 w )Z» + eA"] [{d,xt) Xs ~ fei) xt] 

+(g/2) [- (P X $) X° S W- + (d» X ° s ) X + S W- + h.c] , (2) 

where, e = y/Aira with a being the fine structure constant, g is the SU(2)^ gauge 
coupling constant, and gz = e/(swCw)- The symbols sw and c\y stand for sin 9w 
and cos9\y, respectively, with Q w being the Weinberg angle. 

2.2 Supersymmetric model 

Supersymmetry (SUSY) is the symmetry that relates particles of one spin to other 
particles that differ by half a unit of spin pj. A new particle called superpartner 
is hence introduced for each SM particle in the SUSY model. It is known that the 
chiral symmetry guarantees the smallness of fermion masses. Since SUSY relates 
fermions to bosons, not only the smallness of fermion masses but also that of scalar 
masses are guaranteed, and the hierarchy problem of the SM disappears. In the 
SUSY model, if the R-parity is conserved, the lightest superpartner (LSP) is a good 
candidate for dark matter. One of the most plausible candidates for the LSP is the 
neutralino (xf) which is a linear combination of superpartners of U(l) and neutral 
SU(2)x, gauge bosons and neutral Higgs bosons. On the other hand, a new charged 
particle (xf) 1S a ^ so predicted, namely, the chargino, which is a linear combination 
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of superpartners of the charged SU(2)l gauge boson and the charged Higgs boson. 
The SUSY model therefore predicts fermionic new particles x°f an d Xf- 

In this model, interactions needed for our simulation study have the form: 

£> = ~9zXfY 1 (N l P l + N R P R ) Xf Z, - 9 XfY 1 (C L P L + C R P R ) X ° F W~ + h.c, (3) 

where, Pl and P R are chirality projection operators. Coefficients Nl, N r , Cl, and 
C R in front of the operators are determined by the mass matrices of neutralinos 
and charginos [5J, which depend on the details of the scenario. The values of the 
coefficients adopted in our simulation study are given in the next section. 



2.3 Littlest Higgs model with T-parity 

The littlest Higgs model with T-parity is based on a non-linear sigma model describ- 
ing SU(5)/SO(5) symmetry breaking [7], and the Higgs boson is regarded as one of 
the pseudo Nambu-Goldstone bosons arising from the breaking. The global symme- 
try SU(5) is not exact and is slightly broken due to the existence of explicit breaking 
terms, which are specially arranged to cancel quadratically divergent corrections to 
the Higgs mass term at 1-loop level. The quadratically divergent corrections ap- 
pear, at most, at 2- loop level and the scale of the new physics can be as large as 
10 TeV without the fine-tuning on the Higgs mass term [10], thereby solving the 
little hierarchy problem [9j. Additionally, the implementation of the Z2 symmetry 
called T-parity to the model has been proposed to evade severe constraints from 
electroweak precision measurements [II]. Due to the discrete symmetry, the lightest 
T-parity odd particle, which is the heavy photon (xv)> * s a g°°d candidate for dark 
matter. On the other hand, the charged new particle (xv) which decays into Xv an d 
W is the heavy W boson. Both Xv an d Xv are massive vector bosons, and acquire 
their masses through the SU(5)/SO(5) braking. 

In this model, interactions relevant to our simulation study are given by 

£ = ig[ (c w Z + s w A)^Xv»{d >M Xv U -d' / Xv' 1 ) 
-{c w Z + s w A),x v , (d'xV - d u xV) 
+d, (c w Z + s w A) u (xv"Xv ~ X7XV) 
+shW; X vu {d^Xv ~ d u Xv") ~ shW+x vu {^Xv ~ ^x7) 

+ SH d tl w+ {x° v u x7 - xtyxv) - *hW^ Xvu (0"x£ - STjff) 

+shW~ X °vu {d'xV ~ d v xV) s B d^W- { X °vxV ~ X%X?) ] . ( 4 ) 
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m x ± [GeV] 


m x o [GeV] 


Cross section [fb] 


V~s [GeV] 


Point I 


232 


44.0 


40 & 200 


500 


Point II 


368 


81.9 


40 & 200 


1000 



Table 2: Representative points used in our simulation study. 



where sh = sin^ with Oh being the mixing angle between neutral heavy gauge 
bosons and determined by the mass matrix of the bosons. The value of Oh used in 
our simulation study is also given in the next section. 



3 Simulation framework 

In this section, we summarize the simulation framework such as representative points 
used in our analysis, strategy to discriminate the new physics models discussed in 
the previous section, and tools used in the simulation study. 



3.1 Representative points 

Mass spectrum of the WIMP dark matter (%°) and the new charged particle (x^) 
used in our analysis is shown in Table [2j This mass spectrum is adopted in all the 
new physics models. Though the three new physics models predict different cross 
section values for P arr production, we use a common value for the cross section 
with 100% branching ratio for the decay x ~~ >* X°W ± . Two cross section values are 
considered, 40 and 200 fb, as shown in Table [2J We therefore call the models the 
inert Higgs-like (IH-like), supersymmetric-like (SUSY-like), and littlest Higgs-like 
(LHT-like) models, respectively, in the following discussions. 

In the inert Higgs model, the structures of interaction vertices x + ~X~-Z{l) an d 
X ± ~X°-W T are completely fixec*]. On the other hand, in the supersymmetric model, 
there are parameters to be fixed for the vertices: N L ( R ) and Cur)- With the SUSY 
parameters, m = 5 (10) TeV, M l = 44.5 (81.0) GeV, M 2 = 234 (369) GeV, fi = 1 (1) 
TeV, and tan/3 = 10 (10) at the TeV scale, the masses of the lightest neutralino (%°) 



1 The cross section for the x production in the inert Higgs model is 3.51 (6.85) fb at Point I (II). 
These cross sections are much smaller than the corresponding cross sections in the supersymmetric 
model and the littlest Higgs model with T-parity (see footnotes 2 and 3). This is because the 
production cross section is proportional to /3 3 , in addition to the fact that the production rates for 
scalar particles are usually smaller than those for spin 1/2 fermions or vector bosons. 
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and chargino (X ± ) turn out to be 44.0 (81.9) GeV and 232 (368) GeV, respectively, 
at Point I (II). Using the parameters, the ratio of the coefficients between Nl and Nr 
is determined to be N L /N R = 0.992 (1.00), while C L /C R is 1.36 (1.31). We adopt 
these coefficients in the SUSY-like models. As in the supersymmetric model, there is 
a parameter to be fixed for the vertices in the littlest Higgs model with T-parity: Oh- 
By choosing the vacuum expectation value of the SU(5)/SO(5) symmetry breaking 
to be 375 (580) GeV, we can adjust the masses of x° an d X ± t° be 44.0 (81.9) GeV 
and 232 (368) GeV. With this vacuum expectation value, the angle Oh is determined 
as tan#H = —0.0525 (—0.0246), and we use these values in the LHT-like modej^l 

3.2 Simulation strategy 

Since the dark matter will escape without detection, the measurement of the new 
physics models at the TeV scale (IH-like, SUSY-like, and LHT-like models) is not 
straightforward. In the paper, in order to discriminate the new physics models, we 
focus on the following three physical quantities, (i) the energy distribution of the 
W boson, (ii) the angular distribution of the new charged particle x ± i an d (hi) the 
threshold behavior of the cross section for the X s1 P arr production. These quantities 
are relevant to kinematics of the process and spin information of the new charged 
particles. In this subsection, we discuss how measurements of these quantities work 
for discrimination of the new physics models. 

3.2.1 Energy distribution of W 

Solving the kinematics of the new physics process e + e~ — > x + X~ ~~ >" X°X°W + W~ , 
we find the maximum and the minimum of the W energy (E mSiX and i? m in) given by 

E min = l x ±E* w - /3 x ±7 x ±p^, (5) 

where (3 X ± (7 X ±) is the [3 (7) factor of x^ 1 m the laboratory frame, while (p^) is 

the energy (momentum) of the W boson in the rest frame of x ± - The energy is 

given as (M^± + — M^ )/(2M x o). As a result, both masses of x^ 1 an d X° can be 

estimated from the edges of the distribution of the reconstructed W boson energy. 

2 Cross section for the x ± production in the supersymmetric model is 414 (201) fb at Point I 
(II). 

3 Cross section for the x ± production in the littlest Higgs model with T-parity is 364 (693) fb at 
Point I (II). 
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3.2.2 Angular distribution of the x ± production 

The production angle of the new charged particle x ± can be reconstructed up to 
two-fold ambiguity from the reconstructed W boson momenta. The distribution of 
the reconstructed x ± production angle allows us to investigate the property of x ± , 
because it depends on the spin of x ± - The angular distribution in each case of the 
new physics models (IH-like, SUSY-like, or LHT-like model) turns out to be 



the X s1 momentum and the beam axis. As demonstrated in the following sections, 
the angular distribution turns out to be a powerful tool to discriminate the new 
physics models. 

3.2.3 Threshold behavior of the x ± production 

Since the x ± P a i r production occurs in energetic e + e~ collision through s-channel 
gauge boson exchanges, the total angular momentum along the beam axis in the 
initial state is one. The orbital angular momentum, therefore, has to be one (P- 
wave) when x ± is a scalar particle, which leads to the behavior of the cross section 
a oc (s — 4M 2 ± ) 3//2 in the threshold region s ~ AM\. On the other hand, when 
the x ± is a Dirac fermion, it can be produced with the S-wave, leading to the 
threshold behavior a oc (s - AMl ± ) l l 2 . In the case of the vector x ± , the situation 
is more complicated. Since the x ± m the littlest Higgs model is a gauge boson, 
the production vertex is coming from gauge self-interactions. In addition, there is 
also a vertex between the SM gauge bosons and would-be Nambu-Goldstone bosons 
absorbed in the longitudinal mode of x ± - 111 both cases, the final state with the 
total spin 1 cannot be composed by the vertices alone, which leads to the threshold 
behavior a oc (s — 4M^±) 3 / 2 in the vector production. The threshold behavior of 
the x + production can therefore be used to discriminate the SUSY-like model from 
the rest. 




(for IH-like), 

(for SUSY-like), (6) 

(for LHT-like), 



where x = s/M; 



r2 ± with s being the center of mass energy and 9 is the angle between 
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3.3 Simulation tools 



3.3.1 Event generation 

For both Points I & II, we generated signal events by using the Physsim [12] pack- 
age. In this package, helicity amplitudes are calculated using the HELAS library [13J, 
which deals with the effect of gauge boson polarizations properly. Phase space inte- 
gration and the generation of parton four- momenta are performed by BASES/SPRING [13] ■ 
Parton showering and hadronization are carried out by using PYTHIA6.4 [15], where 
final-state tau leptons are decayed by TAUOLA [16J in order to handle their polar- 
izations correctly. 

For Point I, we generated SM background events by using the matrix element 
generator WHIZARD [TF] 1.40 with PYTHIA 6.205 for the hadronization. It is the 
standard sample for Letter of Intent (Lol) [TSJ [Tj| study of ILC detector concepts, 
covering the whole SM processes with 12 million events in total. In contrast, for 
Point II, SM background events were generated using the Physsim package. 

In all the generated samples, initial-state radiation and beamstrahlung effects 
are included. We ignore the finite crossing angle between the electron and positron 
beams and assume no initial beam polarizationdi. 



3.3.2 Detector simulation 

For Point I, a full simulation code [201 [21], developed for the International Large 
Detector (ILD) [18], is used for the Monte-Carlo (MC) simulation and event re- 
construction. The standard geometry for the ILD Lol study is used for the detector 
simulation. The geometry includes a time projection chamber with silicon devices for 
tracking and vertexing, and highly granular electromagnetic and hadronic calorime- 
ters for particle flow calorimetry along with a 3.5 Tesla magnetic field. The central 
part of the reconstruction is a particle flow algorithm [22], which reconstructs in- 
dividual charged and neutral particles from tracks and calorimeter clusters. The 

4 In general signal and background cross sections depend on beam polarization combination. In 
this study, however, we use no beam polarization so as to keep our study as model-independent 
as possible. If some enhancement is observed for a certain beam polarization combination, we can 
certainly use it to increase our signal statistics. In the case of the maximum enhancement, where 
the signal process is through a single e~ and e + polarization combination, the enhancement is a 
factor of 2.26 for the nominal beam polarizations of 80% in electrons and 30% in positrons. Since 
most of the background processes are enhanced with the left-handed electrons, the background can 
be significantly suppressed if the signal process favors the right-handed electrons. 



9 



Detector 



Performance 



Coverage 



Vertex detector 5 b < 5 © 10/pfi sin 3/2 9 (/an) | cos 6\ < 0.93 

Central drift chamber Sp t /p^ < 5 x 10~ 5 (GeV/c) -1 | cos9\ < 0.98 

EM calorimeter cr E /E = 17%/ \[E © 1% | cosfl| < 0.99 

Hadron calorimeter a E /E = 45%/ y/E © 2% | cos6>| < 0.99 

Table 3: Detector parameters used in the Point II study. 



reconstructed particles are clustered into 4-jet configuration using the Durham al- 
gorithm [23]. A neural-net based flavor tagging algorithm [21] is applied to the jets 
after the jet clustering. 

For Point II, we use a fast simulator code [25], which implements the GLD ge- 
ometry and other detector performance related parameters [26]. In the simulator, 
hits by charged particles at the vertex detector and track parameters at the central 
tracker are smeared according to their position resolutions, taking into account corre- 
lations due to off-diagonal elements in the error matrix. Since calorimeter signals are 
simulated in individual segments, a realistic simulation of cluster overlap is possible. 
Track-cluster matching is performed for the hit clusters in the calorimeter in order 
to achieve the best energy flow measurements. The resulting detector performance 
in our simulation study is summarized in Table [3j 



4 Results from simulation study 

In this section, we present results from our simulation study for e + e~ — > x + X~ ~~ ^ 
X°X°W + W~ process in the case of the IH-like, the SUSY-like, and the LHT-like 
models. We take two cross section points: a s = 200 fb and 40 fb as examples. The 
simulation was performed at y/s = 500 GeV for Point I, and 1 TeV for Point II. An 
integrated luminosity £ int = 500 fb -1 is assumed in each point for all the following 
study except for threshold scans. 

4.1 Study for Point I with yfs = 500 GeV full simulation 
4.1.1 Signal Selection 

Point I employs m x ± = 232 GeV and m x o = 44.0 GeV, which can be investigated 
at the y/s = 500 GeV ILC. We select signal events with both Ws decaying into two 
quarks (qqqq events), whose branching fraction is about 46%, since the W energies 



10 



must be fully reconstructed for the mass determination and the production angle 
reconstruction. The target event topology is thus 4-jets with missing momentum. 
All SM processes with up to 6 particles in the final state are used in the analysis 
as background. The dominant contribution is the W^-pair production with fully 
hadronic decays, the WW Z processes with the Z decaying to a neutrino pair, the 
top-pair production with one W decaying leptonically, and 77 —> WW processes. 
The SM Higgs (ZH, vvH) processes with mjj = 120 GeV and semi-leptonic signal 
processes are also included. 

To reject a major part of the SM and the semi-leptonic decay background, we 
applied primary selection cuts to all samples as follows, (i) The number of tracks 
should be larger than 20 and each jet has to contain at least two tracks in order to 
eliminate pure leptonic events, (ii) The visible energy of the event, E v - ls , should be 
between 80 and 400 GeV, which can remove most of 2-photon and 2, 4, and 6-quark 
events, (iii) Each jet should have a reconstructed energy of at least 5 GeV and a polar 
angle 9 fulfilling | cos^l < 0.99 to ensure proper jet reconstruction, (iv) The distance 
parameter of the Durham jet algorithm [23] for which the event changes from 4-jet 
to 3-jet configuration, 1/34, should be larger than 0.001 in order to reject most of 2-jet 
events, (v) No lepton candidate with an energy larger than 25 GeV is allowed in 
order to suppress semi-leptonic events, (vi) | cos^l of the missing momentum should 
be smaller than 0.9 and | cos#| summed up for all jets should be smaller than 2.6 
in order to eliminate most of the SM events which are concentrated in the forward 
region, (vii) The neural-net output of 6-tag probability summed up for all jets should 
be smaller than 1 to remove events with h quarks. 

After the primary selection, a constrained kinematic fit [27J, which requires the 
two dijet masses of the event to be equal, was performed on each event. All three pos- 
sible jet pairings are tested and the pairing with the least \ 2 value for the kinematic 
fit is selected for the following analysis. 

Secondary selection cuts were applied after the kinematic fit as follows, (viii) 
The kinematic fit constraining the two dijet masses to be equal should converge for 
at least one jet pairing to ensure integrity of the fit result, (ix) The di-jet mass 
obtained by the kinematic fit should be between 65 and 95 GeV to select two-W 
events. 

The effect of these cuts is summarized in Table H] and the distributions of some 
cut variables are shown in Figure [TJ Clear peaks at di-jet masses of 80 GeV can 
be seen in the signal distributions of Figured] (ix), which are from two W bosons. 
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Process 


# of events 


# of events after cuts 




IH-like (hadronic decay) 


46,815 


27,837 


Signal 


SUSY-like (hadronic decay) 


45,550 


26,578 




LHT-like (hadronic decay) 


46,644 


27,631 




IH-like (other decay) 


53,186 


122 


Model bkg. 


SUSY-like (other decay) 


54,462 


104 




LHT-like (other decay) 


53,355 


212 




qqqq (WW, ZZ) 


1.88 x 10 6 


3,218 




qqiv (WW) 


2.35 x 10 6 


1,883 




qqqqvv (WW Z) 


4,158 


681 


SM bkg. 


qqqqtv (tt) 


125,205 


626 




77 ->■ qqqq 


26,356 


509 




qq 


6.29 x 10 6 


373 




SM Higgs (120 GeV) 


56,967 


61 




Other background 


3.44 x 10 9 


338 



Table 4: Event numbers before and after the selection cuts, normalized to £ int = 500 
fb _1 and a s = 200 fb in the Point I study. 
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Figure 1: Cut plots with a s = 200 fb, £ int = 500 fb _1 in the Point I study. The 

labels (ii), (vii), (ix) correspond to the cuts described in the text with the same 

labels. Grayed regions are cut out with the selection. 
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Physics model 


a s = 200 fb 


a s = 40 fb 


M x ± (GeV) 


IH-like 

QTTQV lilro 
jUji -llKe 

LHT-like 


232.9 ± 0.1 

ZOZ. 1 I U.l 

232.1 ± 0.1 


231.8 ± 0.4 

OQO O _|_ n ^ 
ZOZ.Z It u.o 

231.5 ± 0.5 


M x o (GeV) 


IH-like 
SUSY-like 
LHT-like 


44.2 ± 0.6 
43.6 ± 0.7 
43.8 ± 0.5 


46.2 ± 1.9 

45.8 ± 2.3 

45.9 ± 1.8 



Table 5: Measurement accuracies for the masses of x an d X° with £ int = 500 fb 1 
in the Point I study. 

Acceptances of signal events after the cuts are 59.5, 58.3, and 59.2 % for the IH-like, 
the SUSY-like and the LHT-like models, respectively. Signal purities after the cuts 
are 77.4, 76.6, and 77.3 % in the a s = 200 fb case and 40.7, 39.5 and 40.5 % in the 
a s = 40 fb case, respectively. 

4.1.2 Mass Determination 

The masses of new particles can be obtained via the energy spectrum of the W boson 
candidates. The energy of the W bosons has upper and lower kinematic limits, from 
which the masses of the new particles can be derived. Figure |2(a) shows the W 
energy spectrum for each model on top of the SM background. Clear edges can be 
seen in the distribution of every model. 

The edge positions are obtained by a fit using an empirical function with kine- 
matical edges. The analysis was done by three steps as follows. (1) Determine shape 
paremeters of the fitting function with a high-statistical sample (about 1 million 
events per model). (2) Determine the edge positions and the normalization factor 
(three free parameters) by a fit to a sample with the signal cross section (200 fb and 
40 fb). (3) Calculate the masses of x ± an d X° with the obtained edge positions. The 
measurement of the edge positions are assumed to be statistically independent. 

Figures E](b) and [2]^c) give the fitting results of a s = 200 fb and 40 fb with 
Ant = 500 fb -1 , respectively. The fitting results are summarized in Table [5j While 
the central values of the fitting results deviate from the expected masses, they can 
be corrected using Monte-Carlo samples in the real experiment. 

In the W energy distributions, we can see a clear difference among three models, 



14 




W energy (GeV) 




50 100 150 200 250 

W energy (GeV) 




W energy (GeV) 



Figure 2: (a)W energy distributions for signal (cr s = 200 fb) and background with 
£j nt = 500 fb -1 in the Point I study. (b),(c) Results of the mass fit for <j s = 200 fb 
and 40 fb after background subtraction, respectively. 
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which may be used for the model separation. However, the difference is considered 
to be coming from the vertex structures of interactions of the specific models and 
not from the general spin structure, thus we do not use this difference so as to keep 
this study model-independent. 

4.1.3 Angular Distribution for x ± Pair Production 

The separation of the three models is possible by comparing the distributions of 
X ± production angles. To derive the production angles, a quadratic equation is 
solved using the masses of new particles and the momenta of the W bosons with 
the assumption of a back-to-back ejection of the x^ pair- The equation gives either 
two solutions which contain one correct production angle or no solutions when the 
discriminant of the equation is negative. The unphysical negative discriminant comes 
from misreconstructing W momenta or imperfect back-to-back condition of the two 
X ± mainly due to initial state radiation. Fractions of 23.9% (IH-like), 20.8% (SUSY- 
like), 23.7% (LHT-like), and 64.4% (SM background) of the events have negative 
discriminant and are discarded before the following analysis. 

Figure E] shows the production angle distributions. One-dimensional results (a)(b) 
show the visible difference among the three models that the IH-like events concen- 
trate in the central region while the SUSY-like and the LHT-like events are almost 
flatly distributed. Two-dimensional results (c)-(f) are actually used to estimate the 
separation power. We compare the two-dimensional production angle distribution 
for one model (dubbed as "dataset") against another model ("template"). Distribu- 
tions (c)-(e) are used as templates for each model after adding the SM background 
(f). Datasets for each model are created by fluctuating each bin of the templates 
with Poisson distribution. To quantify the difference between a dataset of the model 
M D and a template of the model M T , we defined the chi-square x 2 (M D , M T ), the 
reduced chi-square x 2 (M£>,Mr) and the separation power P(Mo, Mr) as 



where Z)j(M) and Tj(M) are the numbers of the dataset and the template events in 
the ith bin of the model M, N = 210 is the number of bins and a(M) is the standard 




X 2 (M D ,M T ) 



P(M D ,M T ) 



X 2 (M D ,M T ) 

N — 1 
X 2 (M D ,M T )-1 
a(M T ) 



(7) 
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Figure 3: Production angle distributions in the Point I study, (a) and (b) show the 
generated and the reconstructed 1-dimensional distributions. Both of two solutions 
of the quadratic equation are included in (b). The difference between (a) and (b) 
reflects the effect of the wrong solution as well as the detector response, (c)-(f) give 
2-dimensional distributions of the IH-like, the SUSY-like, the LHT-like models, and 
the SM background, normalized to a s = 200 fb and £ int = 500 fb _1 . 



17 





T\/T \ H/T 
1V1 D \1V1T 


ttt 


oUoi -like 


T HT lil^ 
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zUU ID 


ttt 




DO 


A Q 




Q.TTQ.V lilro 


oo 




A Q 




LHT-like 


35 


4.9 




40 fb 


IH-like 




8.9 


6.0 




SUSY-like 


7.5 




0.7 




LHT-like 


4.9 


0.8 





Table 6: Expectation value of separation power P between the three models with 
the 2- dimensional production angle distribution with £ int = 500 fb^ 1 in the Point I 
study. 

deviation of the x 2 {M, M). Since we use a high-statistics sample (1 million events for 
each model) for the template, the effect of the MC statistics of the templates can be 
ignored. The template distributions are normalized to the integral of the data events 
before calculating the x 2 value. Figure H] shows the obtained x 2 distribution with 
10,000 datasets for every combination of the three models. Separation is possible for 
every model with a s = 200 fb, while in the a s = 40 fb case clear separation between 
the SUSY-like and the LHT-like models is impossible. 

Table [6] tabulates the expected values of obtained separation power P. Despite 
the similar angular distribution of the SUSY-like and the LHT-like models, all the 
three models can be identified with a s = 200 fb. In the a s = 40 fb case, the SUSY- 
like and the LHT-like models cannot be separated while the IH-like model can still 
be separated from the other two. These values do not include the effect of the mass 
uncertainty of new particles, which is not significant with < 5% mass uncertainty 
obtained in our mass determination analysis (see Table El). 

4.1.4 Threshold Scan 

Another strategy to distinguish the models is the threshold scan. Figure [5] (a) shows 
how the cross section of each model depends on y/s. A clear difference can be seen 
between the SUSY-like model whose production cross section has the (s — So) 1 ^ 2 
dependence and the other two whose cross sections have the (s — s ) 3 ^ 2 dependence 
where so is the threshold energy, which is twice the mass of the x ± - 

To estimate the separation power with the threshold scan, we performed a toy- 
MC study, in which the measured cross section was fluctuated using the expected 
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Figure 4: x 2 (Md, Mt) distributions with £ int = 500 fb 1 in the Point I study. 
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Figure 5: (a) Dependence of the cross section on the center of mass energy, normal- 
ized to a s = 40 fb in the Point I study. Error bars are given assuming £j nt = 50 
fb -1 data at each point. (b),(c) Results of the \ 2 fits for the (s — Sq) 1 ' 2 case and the 
(s — so) 3//2 case, respectively. 

signal and background statistics obtained from the full-MC study. The cut efficiency 
and the background cross section are assumed to be identical to the 500 GeV case 
for any yfs . We performed a scan of three points: \fs = 470, 485, and 500 GeV, 
each with £- mt = 50 fb -1 . The cross section is scaled so that all the three models 
have a s = 40 fb at 500 GeV. 

For the separation, we calculate the \ 2 value of the fit of 

a(s, n) = a(s - s ) n , n = 1/2, 3/2 (8) 

where a and So are the free parameters for each model. Figures 0(b) and (c) show the 
X 2 distributions. With the n — 1/2 fit (b), good separation is obtained between the 
SUSY-like and the other two models. For example, 92.0% of the SUSY-like events 
are within x 2 < 3 while 5.7 and 2.1% of the IH-like and LHT-like events remain in 
the same \ 2 region. The n = 3/2 fit (c) does not have significant separation power. 
Separation between the IH-like and the LHT-like models is almost impossible by the 
threshold scan. 

Since the SUSY-like model can be separated from the LHT-like model with the 
threshold scan and the IH-like model can be separated from the SUSY-like model 
with the production angle distribution, the three models can be separated from each 
other with combining the two methods even in the a s = 40 fb case. 
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4.2 Study for Point II with yfs = 1 TeV fast simulation 

Since most of the analysis procedure is the same as in the Point I study, we mainly 
focus on the difference and the result of the Point II study in this subsection. 

4.2.1 Signal selection 

Point II (m x ± = 368 GeV, m x o = 81.9 GeV) is not accessible with yfs = 500 GeV 
ILC, so we use 1 TeV fast simulation for the Point II study. As in the Point I study, 
hadronic decay modes of W bosons have been used to select the signal process. All 
events were reconstructed as 4-jet events by adjusting the cut on y-values. In order 
to identify the two W bosons from x ± decays, two jet-pairs were selected so as to 
minimize a x 2 function, 

X 2 = r c M wl - tv M w f/a 2 Mw + r c M W 2 - tT M w f/a 2 Mw , (9) 

where iec Mwi(2) is the invariant mass of the first (second) 2-jet system paired as a W 
candidate, tr M w is the true W mass (80.4 GeV), and <Jm w is the resolution for the 
W mass (4 GeV). We required \ 2 < 26 to obtain well-reconstructed events. Since 
X°'s escape from detection resulting in missing momentum, the missing transverse 
momentum ( miss p<r) of the signal peaks at around 175 GeV. We have thus selected 
events with miss p T above 84 GeV. The numbers of events after the selection cuts 
are summarized in Table [7J Leptonic decay in x ± P arr production and SM Higgs 
backgrounds are not included in the Point II study. These backgrounds are expected 
to be small according to the Point I study. 

4.2.2 Mass determination 

Procedure of the mass determination is almost the same as in the Point I study. The 
masses of x° an d X s1 were determined from the edges of the W energy distribution 
shown in Fig. [51 After subtracting the backgrounds, the distribution was fitted with 
a line shape determined by a high statistics signal sample. The fitted masses of x° 
and x ± with £ mt = 500 fb _1 are summarized in Table [HI The masses of x ± an d X° 
are obtained with accuracies of better than 0.3% and 1.5%, respectively, for a s = 200 
fb. For a s = 40 fb, the measurement accuracies of x ± an d X° are 0.5-1% and 3-6%, 
respectively. 
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Table 7: The number of events before and after the selection cuts, normalized to 
dm = 500 fb" 1 and a s = 200 fb in the Point II study. 





Physics model 


o s = 200 fb 


a s = 40 fb 


M x ± (GeV) 


IH-like 
SUSY-like 
LHT-like 


367.4 ± 0.9 

368.5 ± 0.8 
367.5 ± 0.6 


366.5 ± 3.4 
370.7 ± 2.8 
367.2 ± 2.0 


M x o (GeV) 


IH-like 
SUSY-like 
LHT-like 


81.2 ± 1.1 
81.6 ± 1.1 
82.1 ± 0.8 


80.5 ± 4.7 
82.5 ± 4.5 
84.0 ± 2.7 



Table 8: Measurement accuracies for the masses of x ± and X° with Ant = 500 fb 
in the Point II study. 
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Figure 6: (a) W energy distributions for the signal (cr s = 200 fb) and background 
with d n t = 500 fb^ 1 in the Point II study. (b),(c) Results of the mass fit for cr s = 200 
fb and 40 fb after background subtraction, respectively. 
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Table 9: Expectation values of separation power P between three models with the 
2- dimensional production angle distribution with £ int = 500 fb -1 in the Point II 
study. 

4.2.3 Angular distribution for pair production 

The model separation was studied using the two-dimensional production angle dis- 
tributions as in the Point I study. Figure [7] shows the one- and two-dimensional 
histograms for the two solutions of the production angle. The angular distributions 
for each physics model were prepared with a high-statistics sample, and normalized 
to C mt = 500 fb -1 . The number of bins in 2-dimensional histograms (N in Eq. (j7]) 
) is 325 (instead of 210 in the Point I study). Figure [8] shows the x 2 distributions 
and Table [9] tabulates the expectation values of separation power P for each physics 
model, which are defined in Eq. (j7j). The physics model can be identified confidently 
by using the P values in both of the a s = 200 fb and 40 fb cases. 

4.2.4 Threshold scan 

Threshold scan was also performed with the same procedure as in the Point I study. 
Figure [9] (a) shows the y/s dependence of the cross section of each model. We 
performed a 3-point toy-MC scan of \fs = 750, 800, and 850 GeV with £ int = 50 
fb -1 at each point. The signal cross section a s was scaled to 40 fb at y/s = 1 TeV. 
The cut efficiency and background cross section are assumed to be the same as those 
at y/s = 1 TeV. The \ 2 distributions of fits to Eq. (JSJ) with n = 1/2 and 3/2, shown 
in Figs. M (b),(c), were obtained by the same methods as in the Point I study. Both 
distributions give good separation between the SUSY-like and the other two models. 
If we assume events with \ 2 < 16 as the SUSY-like events for the fitting with the 
power of 1/2, the IH-like and the LHT-like events were disfavored with probability 
of 99.4% and 90.7%, respectively. Here, 92.0% of the SUSY events were selected as 
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Figure 7: Production angle distributions in the Point II study, (a) and (b) show 
the generated and the reconstructed 1-dimensional distributions, (a) shows the true 
distribution, while (b) includes both of the two solutions, (c)-(f) give 2- dimensional 
distributions of the IH-like, the SUSY-like, the LHT-like, and the SM background, 
normalized to a s = 200 fb and C mt = 500 fb _1 . 
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Figure 9: (a) Dependence of the cross section on the center of mass energy, normal- 
ized to cr s = 40 fb at 1 TeV in the Point II study. Error bars are given assuming 
Ant = 50 ftr 1 data at each point. (b),(c) Result of the x 2 fit for the (s — Sq) 1 ^ 2 case 
and the (s — s ) 3 ^ 2 case, respectively. 

the SUSY-like events. On the other hand, if we assume events with x 2 > 13 as the 
SUSY-like events for fit with the power of 3/2, the IH-like and the LHT-like events 
were disfavored with probability of 97.5% and 90.8%, respectively, and 92.3% of the 
SUSY-like events were selected as the SUSY-like events. 

5 Summary 

The WIMP dark matter is one of important candidates predicted in many new 
physics models at the TeV scale, which will be detected at the ILC Interestingly, 
various new physics models predict the existence of the process e + e~ — >■ x + X~ ~~ ^ 
W + W~x°X° i which allows us to measure properties of the dark matter (%°) and the 
new charged particle ix^) with good accuracy. With the use of the process, it is also 
possible to discriminate the new physics models in a model-independent way. We 
have shown that the masses of x° an d X > the angular distribution of x > an d the 
threshold behavior of the x production cross section can be accurately measured at 
the ILC. In fact, it was shown quantitatively that these measurements can be used 
to discriminate the new physics models: IH-like, SUSY-like, and LHT-like models. 

In the study of the benchmark point I, it turns out that the masses of x° an d X 
are determined with accuracies of 5% and 0.2% when the production cross section 
of x^ 1 is a = 40 fb, and 2% and 0.04% when a = 200 fb. The measurement of 
the angular distribution of X s1 enables us to discriminate the IH-like model from 
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the other models, while the SUSY-like model can be discriminated by the threshold 
scan of the process when a = 40 fb. When a = 200 fb, all the new physics models 
can be separated from each other by using only the measurement of the angular 
distributions. On the other hand, in the study of the benchmark point II, the masses 
of x° an d X ± are determined with accuracies of 5% and 0.8% when the production 
cross section of x is er = 40 fb, and 2% and 0.2% when a = 200 fb. The new physics 
models can be discriminated by using the angular distribution even if a = 40 pb. 

In this article, we have shown that new physics models (IH-, SUSY-, and LHT-like 
models) can be discriminated at the ILC On the other hand, it is also true that we 
need to extend the method developed in this article in order to establish a strategy 
for the discrimination in a completely model-independent way. For example, the 
angular distribution of the x + X~ production would be changed if there is a diagram 
in which a new particle (such as selectron in MSSM or heavy electron in LHT models) 
propagates in t-channel. Even if the mass of such a new particle is as heavy as 1 TeV, 
its effect can be sizable in general and the resultant production angle distribution may 
become significantly asymmetric. In addition, if we allow a more generic (Lorentz) 
structure for the x + X Z vertex, the angler distribution may also be affected. In 
these situations, the identification of the W charge [24] becomes very important to 
reconstruct the asymmetric distribution. Moreover, the beam polarization and the 
measurements of the W polarization and the W energy distribution may also play 
an essential role to extract the information on the vertices involving new particles; 
these will help us not only to discriminate new physics models but also to determine 
the properties of the WIMP dark matter in detail. 
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